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Abstract 


This study explores the influence of land altitude on soil organic carbon (SOC) storage in the forest 
soils of Talesh County, Gilan Province. Soil samples were collected from four distinct altitudinal 
ranges (500-1000, 1000-1500, 1500-2000, and 2000-2500 meters above sea level) with comparable 
intrinsic characteristics. Composite soil samples were taken from a depth of 0—35 cm, and the 
concentration of soil organic carbon, combined with soil bulk density from undisturbed samples, was 
used to calculate SOC stocks. 

Statistical comparisons were conducted using a randomized complete block design (RCBD). The Roth 
C model was employed to estimate the organic carbon content of the soils. Climatic data and soil 
organic carbon content at each altitude were utilized to calibrate and validate the model. To assess 
potential variations in SOC storage under changing conditions, diverse climate change scenarios were 
simulated using the Roth C model. 

The findings revealed that as land elevation increased, SOC storage also increased, primarily due to 
higher precipitation and lower temperatures at higher altitudes. The Roth C model demonstrated high 
precision in estimating SOC levels, with correlation and determination coefficients of 0.97 and 0.95, 
respectively. 

The model further suggested that under current climate conditions, soil organic carbon remains stable 
over time. At equilibrium, the carbon input into the soil through plant residues equals the carbon output 
from the soil as carbon gas emissions. 


Keywords:Soil Carbon (SOC), AltitudeGradient, Forest Ecosystems, Carbon Sequestration, 
Roth C model, Climate change Scenarios. 
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حکیده 


این پژوهش به بررسی تأثیر ارتفاع و تغییر اقلیم بر ذخیره کربن آلی خاک‌های جنگلی شهرستان تالش در استان گیلان 
پرداخته است. نمونه‌های مرکب خاک از عمق ۰ - YO‏ سانتی‌متری و از چهار محدوده ارتفاعی (۸۵۰۰-۱۰۰۰ -۱۰۰۰ 
۰ ۱۵۰۰-۲۰۰۰ و ۲۰۰۰-۲۵۰۰ متر بالاتر از سطح دریا) جمع‌آوری و غلظت کربن ui‏ و ple‏ ویژگی‌های خاک 
در آن‌ها اندازه‌گیری گردید. نتایج OUS‏ می‌دهد با افزایش ارتفاع از سطح دریاء ذخیره کربن ST‏ خاک به دلیل افزايش 
بارندگی و کاهش دماء افزایش می‌یابد. بیشترین ذخیرہ مربوط به ارتفاعات ۲۰۰۰ تا ۲۵۰۰ متر (۴۶۔۹۷ تن در هکتار) و 
کم‌ترین آن مربوط به ارتفاعات ۵۰۰ تا ۱۰۰۰ متر (۴۴.۲۳ تن در هکتار) بود. مدل روتامستد Roth C)‏ به‌طور دقیق و 
با ضریب همبستگی و تبیین بالای ۰.۹۵ ذخیره کربن آلی خاک‌ها را تخمین زد. در شرایط فعلی. کربن ورودی و 
خروجی خاک در تعادل و تقریباً معادل ۱.۲۹ تا ۲۷۶ تن در هکتار در سال است. طبق سناریوهای مدل‌سازی, کاهش 
بارندگی (۲.۱۵ میلی‌متر در هر ۱۰ سال) و افزایش دما (۰.۴سانتی گراد)» Gel‏ کاهش ۳۳.۵۴ تا ۴۰.۰۱ درصدی ذخیره 
کربن آلی در همه طبقات ارتفاعی. به‌ویژه در نقاط مرتفع‌تر خواهد شد. در این حالت. خاک به منبعی برای تولید گاز 
کربنیک تبدیل می‌شود. برعکس, افزایش بارندگی و کاهش دما باعث افزایش ۱۳۸۴ تا ۱۸.۹۴ درصدی ذخبره کربن 


آلی در همه طبقات ارتفاعی» به‌ویژه در نقاط مرتفع‌تر خواهد شد. در این حالت» خاک به مخزنی برای نگهداری گاز 
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کربنیک اتمسفری به شکل هوموس تبدیل می‌شود. نتایج OUS‏ می‌دهد مدل‌های شبیه‌سازی دینامیک کربن آلی SLE‏ 
ابزارهای دقیقی برای پیش‌بینی و بررسی تأثیر تغییر اقلیم بر این ذخاین به‌ویژه در خاک‌های جنگلی هستند. این 
پژوهش ast‏ می کند bir‏ ارتفاعات جنگلی و مدیریت اقلیمی برای جلوگیری از کاهش ذخیره کربن و مقابله با 
تغییرات اقلیمی» اهمیت ویژه‌ای دارد. 


کلیدواژەھا:کربن آلی «SU‏ شیب ارتفاعی» اکوسیستم‌های جنگلی»جذب کربنءمدل Roth C‏ سناریوهای ou‏ اقلیم. 


۱-مقدمه 

ماده آلی خاک (SOM)‏ یکی از پیچیده‌ترین اجزاء خاک به‌حساب می‌آید. زیرا از بقایای گیاهی. میکروبی 
و جانوری که در مراحل مختلف تجزیه هستند» تشکیل شده است و این مخلوط ناهمگن با اجزای غیرآلی 
مرتبط است )1992 (Christensen,‏ ماده آلی خاک یک شاخحص مهم در کیفیت خاک است و کمیت و کیفیت 
آن در حفظ حاصلخیزی Coda SY gb‏ خاک سهم بسزایی دارد ( Zhao, Liu, Wang & Zhao, 2015; Ramesh,‏ 
(Manjaiah, Tomar & Ngachan, 2013‏ 
کربن آلی (SOC) SE‏ یک عامل تعیین‌کنندہ در عملکرد اکوسیستم‌های کشاورزی است و مقدار آن در یک 
خاک معین با تفاوت بین میزان کربن ورودی به خاک (بقایای گیاهی و ریشه) و مقدار کربن خروجی از Ol‏ 
(انتشار CO2‏ به هوا) مشخص می‌شود. چندین عامل کنترل‌کننده مانند توپوگرافی ) Kaveh, Mahdian,‏ 
(Parvizi, Sokouti Oskouei & Masihabadi, 2014; Cardinael et al., 2017‏ « آبو هوا ) Mansuri, Karimi,‏ 
(Emamy & Parvizi, 2017; Mufioz-Rojas, Abd-Elmabod, Zavala, De la Rosa & Jordan, 2017‏ « نوع 
خاک )2016 «(Zhao et al., 2015; Kashi Zenouzi, Shafiee & Jafari,‏ عمق نمونه‌برداری )2017 «Li et al.,‏ 
ترکیب کانی‌شناسی )2017 «(Dwivedi et al.,‏ نوع پوشش گیاهی ) Komarov et al, 2017; Rousta,‏ 
(Soleimanpour, Enayati & Pakparvar, 2‏ نحوه استفاده از زمین و شیوه‌های مدیریتی ( Wang et al,‏ 
2015( و برهم‌کنش بین آن‌ها OG‏ کربن ورودی و خروجی را JES‏ و بر مقدار کربن آلی خاک تأثیر 
می گذارند۔ درک عواملی که بر ذخیره‌سازی کربن آلی در خاک SE‏ می گذارند برای مدیریت پایدار در 


1 Soil Organic Matter (SOM) 
2 Soil Organic Carbon (SOC) 
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کاربری اراضی بسیار مهم است. ارتفاع زمین از سطح bos‏ می‌تواند ذخیره کربن آلی را از طریق تغییر 
پارامترهای آب‌وهوایی, تغییر میزان فرسایش خاک» تغییر پوشش ALS‏ و تغییر فرایندهای خاک‌سازی تحت 
تأثیر قرار دهد. بنابراین» ارزیابی تأثیر ارتفاع زمین بر ذخیره‌سازی کربن ull‏ برای برنامه‌ریزی و مدیریت 
پایدار در کاربری اراضی ضروری است. 

ذخیره‌سازی کربن. فرایندی است که در آن دی اکسید کربن جو در مخازن طبیعی یا مصنوعی مانند جنگل‌هاء 
اقیانوس‌ها و خاک ترسیب و ذخیره می‌شود )2008 (Gibbins & Chalmers,‏ جنگل‌ها یکی از بزرگ‌ترین 
مخازن طبیعی کربن در زمین هستند و حدود ۳۰ درصد از ذخیره کربن جهانی را تشکیل می‌دهند. اهمیت 
ذخیره‌سازی کربن در جنگل‌ها به دلیل نقش آن‌ها در کاهش تغییرات اقلیمی در سال‌های اخیر موردتوجه زياد 
قرار گرفته است. به‌علاوه جنگل‌ها خدمات اکوسیستمی متنوعی را ارائه می‌دهند که شامل حفظ تنوع زیستی, 
ذخیره‌سازی آب و حفاظت از خاک است. همچنین با ایجاد معاش برای جوامع محلی. حمایت از گردشگری 
و تولید محصولات چوبی و غیر چوبی به توسعه پایدار نیز کمک می‌کنند. 

مطالعات متعدد نشان داده است که ارتفاع زمین بر توزیع و میزان ذخیره‌سازی کربن آلی در اکوسیستم‌های 
جنگلی "ur‏ می گذارد )2017 (Zhu et al., 2010: Bangroo, Najar & Rasool,‏ شیب و ارتفاع می‌تواند بر 
رطوبت. دما و در دسترس بودن مواد مغذی خاک تاثیر بگذارد که این عوامل به‌نوبه خود بر تجزیه مواد آلی 
و تشکیل کربن آلی در خاک تأثیرگذار است. چندین مطالعه نشان داده‌اند که ذخیره کربن آلی خاک با افزایش 
ارتفاع به دلیل تجمع مواد آلی در لایه‌های بالایی خاک و کاهش سرعت تجزیه افزایش Zhang et al.,) Jub p»‏ 
9 آن‌ها اعلام نمودند که جنگل‌هایی که در ارتفاعات بیش‌تر واقع شدەاند معمولاً دارای دمای کمتر و 
بارندگی بیشتری هستند که باعث رشد بیشتر گیاهان از یک طرف و تجزیه کم‌تر مواد آلی خاک از طرف 
So‏ می‌شود. درنتیجه دی‌اکسید کربن بیشتری از جو جذب و در زیست‌توده گیاه ذخیره می‌شود که بعدا 
به‌صورت بقایای گیاهی به خاک اضافه می‌شود. 

بااین‌حال, مطالعات دیگر نشان داده‌اند که ذخیره کربن آلی خاک با افزایش شیب و ارتفاع به دلیل افزایش 
فرسایش و شستشوی خاک کاهش می یابہد )2015 .(Wang et al.,‏ آن‌ها در مطالعه خود نشان دادند که محتوای 


کربن آلی خاک با افزایش ارتفاع به دلیل دماهای کم‌تر و فصول رشد کوتاه‌تر که منجر به رشد کم‌تر گیاهان 
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می‌شود. کاهش می‌یابد. این نتایج متناقض OUS‏ می‌دهد رابطه بین ارتفاع زمین و ذخیره‌سازی کربن آلی خاک 
پیچیده و بسته به شرایط خاص اکوسیستم جنگلی؛ نوع SE‏ و نوع پوشش‌های گیاهی متفاوت است. 

ارتفاع زمین از سطح dys‏ یک عامل GUIS‏ در تشکیل خاک است و باعث ایجاد یک رابطه قوی ہی کربن 
آلی SE‏ و ویژگی‌های زمین در مقیاس‌های مزرعه‌ای می‌شود )2007 (Papiernik et al,‏ اگرچه چندین 
مطالعه سعی کرده‌اند به رابطه بین کربن آلی خاک و ارتفاع زمین بپردازنده ولی اطلاعات کمی در مورد 
تعاملات led!‏ وجود دارد )2009 .Senthilkumar, Kravchenko & Robertson,‏ دو دلیل برای مقدار بیش ‌تر 
مواد آلی در ارتفاعات کمتر گزارش شدہ است یکی Sy!‏ در ارتفاعات کمتر در مقایسه با ارتفاعات بیش تر 
رطوبت نسبی هوا بیش تر است. دوم این که در ارتفاعات بیش‌تر مواد آلی از طریق رواناب و فرسایش به JUS‏ 
پایین تر منتقل می شود )2002 (Quideau,‏ همچنین سطوح بالاتری از ماده آلی در شیب‌های شمالی در مقایسه 
با شیب‌های جنوبی به دلیل وجود دمای SeS‏ وجود دارد. کربن آلی خاک در نواحی تپه‌ای» عمدتاً تحت تأثیر 
نوع پوشش گیاهی و همچنین ارتفاع قرار می‌گیرد» زیرا ارتفاع تا حد زیادی بر عوامل اقلیمی» عمدتاً دما و 
رطوبت تأثیر می گذارد )2015 (Ahmad Dar & Somaiah,‏ به‌طورکلی» مقدار کربن all‏ خاک با افزایش ارتفاع 
ES‏ کند ) Sinoga, Pariente, Diaz & Murillo, 2012; Choudhury, Fiyaz, Mohapatra & Ngachan,‏ 
6 زرا متغیرهای آب‌وهوایی L‏ ارتفاع تغییر می‌کنند. گاهی اوقات با افزایش ارتفاع» میزان ہارندگی 
سالانه افزايش می‌یابد که این امر فرآیندها؛ ویژگی‌ها و توسعه خاک را کنترل می کند ( Dahlgren, Boettinger,‏ 
(Huntington & Amundson, 1997‏ و AJ i‏ زیست‌توده گیاهی و کربن all‏ خاک را افزايش می دھد ) Sinoga‏ 
al., 2012‏ ). تغییرات ارتفاعی همچنین دما را تغییر می‌دهد که سرعت تجزیه و تجمع کربن آلی خاک b‏ 
کنترل می کند )2016 (Choudhury et al.,‏ گروهی از محققین بیان کردند که در هند با افزایش میزان بارندگی 
میزان کرہن آلی خاک افزایش یافت اما ہا افزایش ches‏ میزان کرہن آلی SE‏ کاهش یافت ) Bhattacharyya et‏ 
GL, 2010‏ بنابراین» اقلیم تأثیر مستقیمی بر نوع و کمیت پوشش گیاهی» سرعت هوازدگی و شدت آب‌شویی 
دارد و درنتیجه کمیت و CAMAS‏ کربن آلی خاک و کیفیت خاک را تعیین می‌کند ) ;1997 Dahlgren et al,‏ 


.(Sinoga et al., 2012 
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مدل‌سازی دینامیک کربن آلی خاک برای پیش‌بینی تغییرات کربن آلی در آینده تحت تأثیر سناریوهای متنوع 
تغیبرات اقلیمی می تواند به درک اثرات متقابل ارتفاع زمین و e)‏ بر ذخیره‌سازی CHS‏ آلی خاک کمک XS‏ 
(Don et al., 2024)‏ مدل‌های کربن آلی به‌عنوان ابزارهای ارزشمندی برای درک تعاملات پیچیدہ بین تغییر 
اقلیم SE‏ و OLS‏ مورداستفاده قرار می گیرند )1987 (Parton, Schimel, Cole & Ojima,‏ ہا درک بهتر این 
فرایندها» توانایی ما در توسعه راهبردهای موثر برای کاهش تغییرات اقلیمی و ترویج توسعه پایدار نیز بهبود 
خواهد یافت )2013 (Aguilera, Lassaletta, Gattinger & Gimeno,‏ 

مدل "Roth C‏ مدل پرکاربردی برای پیش بینی دینامیک کربن آلی خاک در اکوسیستم‌های جنگلی در شرایط 
آب‌وهوایی متفاوت است. این مدل دینامیک کوتاه‌مدت را شبیه‌سازی می کند و ورودی‌های آن ماده آلی SLE‏ 
نرخ تجزیه و فرآیندهای کانی‌ساز را شامل می شود و می تواند برای شبیه‌سازی تغییرات کربن آلی خاک در 
شرایط محیطی و اقلیمی مختلف استفاده شود. 

مطالعات متعددی از این مدل برای بررسی اثرات ارتفاع بر ذخیرەسازی al OUS‏ خاک استفاده کرده‌اند. 
به‌عنوان‌مثال. گروهی از محققین» از مدل Roth C‏ برای شبیه‌سازی دینامیک کربن Ui‏ خاک در یک 
اکوسیستم جنگلی کوھستانی در جنوب غربی چین استفاده کردند و دریافتند که ذخیرہ کربن آلی خاک با 
افزایش گرادیان شیب به دلیل افزایش فرسایش کاهش می‌یابد )2017 Li etal.,‏ برخی دیگر طی مطالعه‌هایی 
در چین با استفاده از مدل Roth C‏ نشان دادند که کربن آلی خاک در طول دهه‌های آیندہ به دلیل تغییر اقلیم 
کاهش می یابد )2011 (Wan, Lin, Xiong & Guo,‏ گروهی مطالعه‌ای را با هدف آزمایش عملکرد مدل Roth‏ 
C‏ با استفاده از داده‌های اندازه‌گیری شده از پنج نوع پوشش جنگلی US ua‏ (شمال ایران) و همچنین 
پیش‌بینی تغییرات کربن آلی خاک تحت سناریوھای مختلف تغییر اقلیم انجام دادند ( Soleimani, Hosseini,‏ 
.Bavani, Jafari & Francaviglia, 7‏ نتایج تجزیه‌و تحلیل داده‌ها OUS‏ داد که مقادیر Sle a S‏ 
شبیه‌سازی شدہ با داده‌های اندازه‌گیری شده همبستگی مناسبی (ضریب تعیین بین ۰/۶۴ تا ۰/۹۱) دارد که 
نشان‌دهنده کارایی مناسب مدل Roth C‏ در این مناطق است. lol‏ اعلام نمودند که تا سال ۲۰۹۹ به‌طورکلی 


ذخایر کربن آلی خاک تحت تمام پوشش‌های موردمطالعه و سناریوهای تغییر اقلیم تعریف شده کاهش یافت. 


' Rothamsted Carbon Model 


۱۵۶ نشریه جغرافیا و مخاطرات محیطی, جلد ۱۳ء شمارہ ۲ پائیز ۱۴۰۳ 


اما میزان کاهش کربن آلی خاک بین مدل‌های آب‌وهوایی. سناریوهای انتشار گاز گلخانه‌ای» دوره‌های زمانی 
و پوشش‌های مختلف گیاهی متفاوت بود. 

اگرچه عامل ارتفاع زمین از سطح دریا نقش مهمی در مقدار کربن آلی خاک‌های جنگل دارد. ولی تأثیر تغییر 
اقلیم بر میزان تغییرات کربن آلی خاک در ارتفاعات مختلف ناشناخته است. لذا هدف تحقیق حاضر 
شہیەسازی رفتار کربن آلی خاک در ارتفاعات مختلف زمین در پاسخ به تغییر اقلیم است. مطالعه تغییرات 
کربن آلی نشان خواهد داد که در آیندہ این اکوسیستمھا چه وضعیتی خواهند داشت» آیا به منبعی برای تولید 
CO2‏ اتمسفری تبدیل می‌شوند و یا مخزنی برای نگهداری دی اکسید کربن BL‏ خواهند ماند. بەطورکلی 
مدل‌ها با ارائه اطلاعات در خصوص نحوه تغیبر کربن آلی ذخیرہ شده در خاک که در پاسخ به عوامل 


۲-مواد و روش‌ها 


۱-۴- منطقه موردمطالعه 


این پژوهش در جنگل‌های شهرستان تالش که در مسیر اسالم به خلخال قرار costo‏ انجام گرفته است (شکل 
-C‏ شهرستان تالش از شهرستان‌های استان OS‏ است که دارای اقلیم معتدل و مرطوب است. میزان بارندگی 
سالانة ol‏ در ایستگاه هواشناسی آستارا (شمالی‌ترین نقطه) ۱۴۴۶ میلی‌متر و در ایستگاه هواشناسی 
رضوان‌شهر (جنوبی‌ترین نقطه) ۱۸۰۰ میلی‌متر است. مختصات جغرافیایی شهرستان تالش ۳۷۰۷۹ درجه‌شمالی 
و ۴۸۹۰ درجه‌شرقی است. جنگل‌های موردمطالعه در ارتفاع ۵۰۰ تا ۳۰۰۰ متری از سطح دریا قرار داشتند و 
مساحت آن‌ها حدود ۱۵۰ هزار هکتار بود. این جنگل‌ها زیستگاه ۱۵۰ گونه درختی (۳۶ گونه بومی) ازجمله 
شمشاد هیرکانی «(Celtis caucasica)‏ لەله خزری (Alnus subcordata)‏ آکاسیای ابریشمی 
.(Platycladus orientalis)‏ بل وط «(Quercus brantii)‏ دودو Juglans nigra)‏ و گردو ) Juglans‏ 


8 می‌باشند. 
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Fig.1. Geographical location of the study area in Iran and Gilan province 


به کمک نقشه‌های جغرافیایی و تحلیل ارتفاعی در محیط نرمافزار QGIS‏ از هر ۵۰۰ متر ارتفاع یک Abb‏ 
ارتفاعی انتخاب و په دلیل نبود ایستگاه هواشناسی مجزا در هر طبقه ارتفاعی برای دریافت داده‌های 
هواشناسی cles)‏ بارش و تبخیر و تعرق)ء از داده‌های شبکه‌بندی شده جهانی با قدرت تفکیک VE‏ در مرکز 
ملی تحقیقات جوی استفاده شد (ترا کلایمت). داده‌های این مجموعه از قدرت تفکیک بسیار بالایی 
پرخوردار بودند که برای این مطالعه بسیار مناسب است. قدرت تفکیک مکانی متغیرهای آب‌وهوایی در این 
پایگاه به‌طور تقریبی Y‏ کیلومتر است. به‌عبارت‌دیگر برای هر چهار کیلومتر در چهار کیلومتر داده‌های 
هواشناسی تخمین زده شده است که OUO‏ از قدرت تفکیک بالا این پایگاه دارد. مقیاس زمانی داده‌های 
استفاده شده در این مطالعه سال‌های ۱۹۵۸ تا ۲۰۲۲ می‌باشد. مطالعه pole‏ چهار طبقه ارتفاعی (۰۰۰-۵۰۰ 


ما ماه YOu Youn‏ متری) را شامل می‌شود. که برای به دست آوردن داده‌های 


pi ۱۵۸‏ 4 جغرافیا و مخاطرات محیطی, جلد ۱۳ء شمارہ ۲ پائیز ۱۴۰۳ 


آب‌وهوای آن‌ها از داده‌های شبکه‌بندی شده جهانی استفاده گردید و شیپ فایل‌های مربوط به هر طبقه 
ارتفاعی به‌صورت مجزا و برحسب مقیاس متر استخراج شد. 

در این پایگاه داده polis‏ دما به‌طور مستقیم وجود نداشت. به همین دلیل به کمک دماهای بیشینه و کمینه در 
آن پایگاه میانگین labs cles‏ واقعی محاسبه گردید (ترا کلایمت). 

ازآنجاکه شیپ فایل‌های اولیه دارای مقیاس متریک بودند و مقیاس داده‌های شبکه‌ی جهانی بر حسب درجه 
LORD eich lli‏ تساه ان po‏ داب - اس ده تل کسی جا تایه از 
als‏ و اسکریپت‌های تعریف شده در محیط برنامه‌نویسی متلب برای هر شیپ‌فایل مرتبط با یک طبقه 
ارتفاعی» مقادیر بارش تبخیر و تعرق و دمای حداکثر و حداقل استخراج و به سری‌های زمانی ماهانه تبدیل 
شدند. درنهایت با گرفتن میانگین سری‌های زمانی ماهانه مقادیر میانگین دماء بارش و تبخیر و تعرق ماهانه 
EM CIT‏ 

مقادیر محاسبه شده میانگین بارش. دما و میزان تبخیر و تعرق ماه‌های مختلف سال به ترتیب در جدول‌های 
٣٢ A‏ مان نان شاه انت uf‏ مس ello Do 8A celal ae Tog‏ سس در لول Loo did‏ 
ماهائه پازاشرھای هراشتاسی برای هر چهار طبقه ارتفاعی از سال ۱۹۵۸ 85ا سال ۲۰۲۲۷ اراته edes‏ است. 


حدول ۱- میانگین بارند گی (میلی‌متر) ماه‌های مختلف سال در طبقات ارتفاعی موردمطالعه. 


Table1- Average rainfall (mm) for different months of the year in the elevation classes under study 


Elev.(m) 


Nov Oct Sep Aug Jul Jun May Apr Mar Feb Jan 

ارتفاع (دی) (بهمن) (اسفند) (فروردین) (اردیبهشت) (خرداد) Gs)‏ (مرداد) (شهریور) (oul) Ge)‏ 
500-0 130.79 147.40 157.41 154.60 64.46 30.52 29.64 23.60 69.85 82.44 187.33 
1000-1500 192.77 173.69 206.96 223.30 112.42 37.78 29.81 28.45 47.57 8393 163.13 
1500-2000 177.37 149.71 188.21 215.93 195.47 172.92 103.29 102.02 122.98 144.21 168.57 
2000-2500 184.74 154.55 202.01 231.22 205.09 175.39 11528 108.62 127.82 157.35 183.20 


جدول ۲- میانگین دما (سانتی گراد) ماه‌های مختلف سال در طبقات ارتفاعی موردمطالعه 


Table2- Average temperature (°C) for different months of the year in the elevation classes under 


study 
Dec Nov Oct Sep Aug Jul Jun May Apr Mar Feb Jan Elev.(m) 
OD — QUD — Gg) (بهمن) (اسفند) (فروردین) (اردیبھشت) (خرداد) (تیر) (مرداد) (شهریور)‎ e») ارتفاع‎ 
17.00 21.90 2827 32.81 36.66 37.19 35.62 4 26.72 21.57 17.36 15.41 500-1000 
11 15.602 22.50 27.46 31.12 31.84 29.1 24.93 19.33 13.28 10.04 5 1000-1500 
7.45 12.52 18.73 23.18 28.87 28.63 26.16 22.49 16.69 10.86 1 5.56 1500-2000 
6.50 11.35 18.35 23.02 26.76 26.80 24.53 20.19 14.56 8.38 5.10 6 2000-2500 
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جدول ۳- میانگین تبخیر تعرق (میلی‌متر) ماه‌های مختلف سال در طبقات ارتفاعی موردمطالعه. 


Table3- Average evapotranspiration (mm) for different months of the year in the elevation 


classes under study 
Nov Oct Sep Aug Jul Jun May Apr Mar Feb Jan Elev.(m) 
(oun Ge) Gere) (مرداد)‎ (s) ارتفاع (دی) (بهمن) (اسفند) (فروردین) (اردیبهشت) (خرداد)‎ 


111 126 1365 114.5 107 101.3 88.14 83.47 82.22 89.25 99.2 500-1000 
106 123 139.3 121.7 106 98.34 85.36 82.22 81.78 84.97 97.1 1000-1500 
878 102 131.1 105.7 98.8 9421 84.37 62.37 56.54 73.42 84.0 1500-2000 
67.2 70.1 83.94 84.84 73.6 59.84 51.92 43.34 39.76 46.47 56.3 2000-2500 


جدول ۴- میانگین بارندگی دما و تبخیر و تعرق سالیانه در طبقات ارتفاعی موردمطالعه. 


Table4- Average annual rainfall, temperature, and evapotranspiration in the elevation classes 
under study 


ارتفاع بارش دما تبخیر تعرق 
Evapotranspiration Temperature Rainfall Elevation‏ 
(mm) (°C) (mm) (m)‏ 
500-1000 1247.26 26.93 1238.38 
1000-1500 148.10 20.45 1231.15 
1500-2000 1937.41 17.28 1055.05 
2000-2500 2052.01 15.83 740.76 


همان‌طور که در جدول‌ها دیدہ می شود با افزایش ارتفاع منطقه متوسط بارندگی سالیانه روند افزایشی و 
میانگین درجه حرارت سالیانه روندی کاهشی دارد. 

۲-۲-نمونه‌برداری خاک و اندازه گیری‌های آزمایشگاهی 

E‏ تأثیر ارتفاع منطقه از سطح دریا بر ذخیره کربن آلی SE‏ چهار طبقه ارتفاعی مختلف در امتداد 
شیب شمالی منطقه موردمطالعه» شامل ارتفاعات ۵۰۰0-۱۰۰۰ ۰٠-۱۵۰۰‏ ۱۵۰۰-۲۰۰۰ و ۲۰۰۰-۲۵۰۰ 
متری از سطح دریا انتخاب شدند. 

در هر طبقه ارتفاعی, حدود ۱۵ هکتار از اراضی جنگلی انتخاب و به سه قسمت پنج هکتاری تقسیم شد و 
هر قسمت به‌عنوان یک تکرار در نظر گرفته شد. سپس, در هر تکرار ۵ پروفیل حفر گردید و در هر پروفیل 
از عمق ۰-۳۵ سانتی‌متری خاک یک نمونه تهیه گردید (به دلیل سنگلاخ بودن خاک‌ها امکان نمونه‌برداری از 
اعماق بیش‌تر میسر نگردید). سپس نمونه‌های خاک به دست آمده از پروفیل‌های مختلف هر تکرار باهم 
ترکیب شدند و یک نمونه مرکب خاک تهیه گردید به‌طوری‌که در هر طبقه ارتفاعی ۳ نمونه مرکب خاک 


وجود داشت. با توجه به وجود چهار طبقه ارتفاعی در این پژوهش تعداد ۱۲ نمونه مرکب خاک جامعه 


۱۴۰۳ جغرافیا و مخاطرات محیطی, جلد ۱۳ء شماره ۲ پائیز‎ 4 pi Me 


آماری این تحقیق را تشکیل می‌دادند. همچنین ۱۲ نمونه دست‌نخورده (عمق ۰-۱۰ سانتی‌متر) نیز برای 
اندازه‌گیری جرم مخصوص ظاهری خاک تهیه گردید. 

پس از اتمام مرحله نمونه‌برداری» نمونه‌های خاک هوا خشک شدند و پس عبور از الک دو میلی‌متری برای 
اندازه‌گیری ویژگی‌های شیمیایی و فیزیکی به آزمایشگاه منتقل شدند. در آزمایشگاه کربن آلی به روش 
واکلی و بلک )1934 (Walkley & Black,‏ بافت خاک به روش هیدرومتر )2002 (Gee & Or,‏ و جرم 
مخصوص ظاهری خاک به روش استوانه )1986 (Blake & Hartge,‏ اندازه گیری شد. 

۳-۲-شبیه‌سازی رفتار کربن آلی خاک در سناریوهای مختلف آب‌وهوایی 

aly‏ مدل‌سازی ahy‏ کرین آلی کاک در ستاریوهای مختلف آب‌وهوایی از YA dead Roth 0 dae‏ استفاده 
ET‏ 

در این مدلء منحصراً به عوامل خاکی توجه نشدہ و مقدار کربن گیاهی که سالیانه از طریق زیتودہ 
گیاهی به خاک اضافه می‌شود. در مدل گنجانده شده است. با استفادہ از این مدل می‌توان چگونگی تجزیه 
بقایای گیاهی وارده به خاک را در شرایط آب‌وهوایی مختلف مشخص نمود و تأثیر al‏ را بر کربن آلی 
خاک تعیین نمود )2004 .(Ponce-Hernandez, Koohafkan & Antoine,‏ 

Roth C ۱-۳-۲-ویژگی‌های مدل‎ 

مدل Roth C‏ برای دینامیک کربن آلی در خاک‌های غیر SUL‏ طراحی شده است و تأثیر نوع خاک des‏ 
رطوبت و پوشش گیاهی را با استفاده از یک دوره ماهانه بر میزان کربن آلی خاک و کربن زیست‌توده 
میکروبی موردبررسی قرار می‌دهد )2017 (Kaczynski, Siebielec, Hanegraaf & Korevaar,‏ در مدل Roth‏ 
C‏ مواد آلی به پنج جز اساسی» شامل ماده آلی خنشی OM)‏ ماده گیاهی قابل‌تجزیه (DPM)‏ ماده 
گیاهی مقاوم به تجزیه (RPM)‏ زیست‌توده میکروبی (BIO)‏ و مواد آلی هوموسی شده (HUM)‏ تقسیم 
می شود )2008 Jenkinson & Coleman,‏ (شکل (Y‏ دو جزء ماده گیاهی قابل تجزیه و ماده گیاهی مقاوم به 
تجزیه مربوط به BLS ChE‏ ورودی به خاک بوده و دو جزء زیست‌توده میکروبی و مواد آلی هوموسی 


شده مربوط به مواد آلی SLE‏ است )2010 (Nemoto,‏ 
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Fig.2. Structure of the Roth C model (Jenkinson & Coleman, 2008)‏ 
برای محاسبه کربن آلی ذخیرہ شده در SLE‏ بر حسب تن بر هکتار از رابطه O)‏ برای هر چهار 


(Nieto, Castro & Fernández-Ondoíio, 2013; Qiu et al., 2013) طبقه ارتفاعی استفاده شد‎ 


CS= 10000 x %SOC x Bd x d :۱ رابطه‎ 


که در این رابطہ:08) کسربن آلی ذخیرہ شدہ در خاک بر حسب کیلوگرم بر هکتار:۹00 درصد 
کربن Bdusk (JE‏ جرم مخصوص ظاهری خاک بر حسب گرم بر سانتی متر مکعب diy‏ عمق نمونه‌برداری 
پر حسب سانتی‌متر است. 

۲-۳-۲ -متغیرهای ورودی مدل Roth C‏ 

بارندگی و تبخیر و تعرق ماهانه: داده‌های بارندگی (میلی‌متر) و تبخیر (میلی‌متر) از سطح آزاد برای 
محاسبه کمبود رطوبت خاک سطحی (TSMD)‏ به‌کار گرفته می‌شوند و دست‌یابی بے آن‌ها برای محاسبه 
کمبود رطوبت خاک سطحی آسان‌تر از اندازه‌گیری ماهانه کمبود رطوبت خاک سطحی واقعی است. اگر 
میزان تبخیر از سطح آزاد در دسترس نباشد. می‌توان متوسط تبخیر و تعرق پتانسیل ماهانه را با 
استفاده از داده‌های هواشناسی تعیین و به کمک رابطه (۲) تبخیر از سطح آزاد را محاسبه و 
شبه‌سازی نم ود )2008 .Jenkinson & Coleman,‏ 


Evaporation = 'Mean potential transpiration'/0.75 ۲ رابطه‎ 


که در این Evaporation akl‏ تبخیر از سطح آزاد و Mean potential transpiration‏ متوسط oz‏ 


d 


pi ۱۶۲‏ 4 جغرافیا و مخاطرات محیطی, جلد ۱۳ء شماره ۲ پائیز ۱۴۰۳ 


متوسط درجه حرارت ماهانه هوا: در این مدل درجه حرارت هوا بیشتر از دمای خاک مورد استفاده قرار 
می گیرد» زیرا درجه حرارت هوا به آسانی قابل دستیابی است. در سطح خاک» متوسط درجه حرارت 
ماهانه هوا نشان‌دهنده متوسط درجه حرارت ماهانه خاک است )2008 Jenkinson & Coleman,‏ 

مقدار رس: درصد رس برای محاسبه مقدار آب در دسترس گیاه که در خاک می‌تواند نگه داشته شود 
به‌کار می‌رود. همچنین» درصد رس بر سرعت تجزیه موادآلی خاک نیز تأثیرگذار است ) & Jenkinson‏ 
Coleman, 2008‏ 

نسبت DPM/RPM‏ نسبت ماده گیاهی قابل تجزیه به ماده گیاهی مقاوم به تجزیه ازجمله شاخص‌های 
تعیین کنندہ نرخ تجزیه مواد گیاهی ورودی به خاک می‌باشد. 

td‏ شا که ای شین los‏ نات دادن سر رعے سفق تا بای میدن غات رو ull‏ که پر اک 
مزروعی سطح خاک دیرتر خشک می شود )2008 Jenkinson & Coleman,‏ 

میزان بقایای گیاهی ورودی به خاک: مقدار بقایای گیاهی که ماهانه به خاک اضافه می‌شود (تن در هکتار 
در ماه) در مدل وارد می‌شود. 

مقدار کود دامی FYM‏ این ورودی متفاوت از ورودی بقایای تازه گیاهی در نظر گرفته شده است 
(تن در هکتار در ماه). 

عمق لایه نمونه‌برداری شده خاک: بەطور معمول خاک سطحی از عمق صفر تا YY‏ سانتیمتر است که 
بسته به شرایط منطقه و نوع پروژه عمق نمونه‌برداری می‌تواند متفاوت باشد. 

مدل Roth C‏ به مقادیر پارامترهای HUM,BIO;RPM,DPM‏ و IOM‏ نیاز دارد که مقدار اولیه 
آن‌ها شناخته شده نیست و فقط مقدار کربن آلی خاک پارامتر قابل دسترس و اندازه‌گیری می‌باشد. به‌منظور 
دستیابی به مقادیر اجزا مختلف کربن آلی خاک. مدل می‌تواند polis‏ این پارامترها را از مقدار کربن آلی 
خاک تخمین و شبیه‌سازی کند. به‌منظور دست‌یابی به نتایج واقعی کالیبره کردن مدل از اهمیت بالایی 
برخوردار است و تجزیه‌های آزمایشگاهی کربن آلی خاک باید دقیق و قابل اعتماد باشد تا به‌توان نتایج 


رضایت‌بخشی را از مدل دریافت نمود. 
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مدل Roth C‏ ویژگی‌هایی مطلوبی مانند قابلیت یادگیری. قابلیت تعمیم ساختار ساده )2010 (Nemoto,‏ 
دسترسین اسان به اطلاعات» سهولت استفاده از مفروضات. تعداد کم پارامترهای مورد SLES‏ 
(Zimmermann, Leifeld & Fuhrer, 2007)‏ و قابلیت اجرا بر ای یک منطقه وسیع ( Shirato & Yokozawa,‏ 
5 را دارد. این ویژگی‌ها کارایی بالایی به مدل, در بعضی از کاربردها مانند سازگاری با شرایط 
منطقه. سطح دقت تخمین‌های مدل و به‌طورکلی در هر جا که نیاز به پیش‌بینی میزان کربن آلی 
شاک با استفاده از بقایای گیاهی و کوهدهی cd‏ می‌بخشد. 
۳-۳-۲-سناریوهای انتخاب شده برای مدل Roth C‏ 
بر اساس پیش‌بینی‌های انجام‌شده. اقلیم کشور 0l ul‏ و همچنین استان گیلان در حال تغییر بوده و در آینده با 
افزایش دما و کاهش بارند گی مواجه خواهد شد )2016 (Kazemi Rad & Mohammadi,‏ آن‌ها با تعریف 
سناریوهای A2‏ و Bl‏ برای مدل‌های گردش عمومی جو MPEHS‏ و HADCM3‏ و با استفاده از مدل 
LARS-WG‏ خروجی این مدل‌ها را ریز مقیاس نموده و میزان تغیبرات دمای هوا و بارندگی را برای هر ۱۰ 
سال پیش‌بینی کردند. براساس نتایج پژوهش daol‏ به ازای هر ۱۰ سال دمای هوا ۰/۴ سانتی گراد افزایش و 
میزان بارش ۲/۱۵ میلی‌متر کاهش می‌یابد. لذا یکی از سناریوهای تعریف شده برای مدل Roth C‏ نتایج 
بەدست stal‏ از مطالعه کاظمی و ھمکاران )2016 (Kazemi Rad & Mohammadi,‏ است که وضعیت ترسیب 
کربن را در شرایط وقوع تغییر اقلیم در منطقه موردمطالعه و ارزیابی قرار می‌دهد. جدول ۵ سناریوهای 
تعریف شده برای مدل Roth C‏ را نشان می‌دهد. 

Roth C سناریوهای در نظر گرفته شده برای مدل‎ -۵ Joa 


Table5- Scenarios considered for the Roth C model 


بازہ زمانی (میلادی) مدیریت سناریو 
Scenario Management Time period‏ 
شرایط کنو 
1958-2022 شرایط کنونی ۱ 
Current conditions‏ 
2022-2102 5 تر "EE‏ 1 
Maintaining the current conditions‏ 
تغییر اقلیم با افزایش ۰/۴ سانتی‌گراد دما و کاهش ۲/۱۵ میلی‌متر 
بارند گی به ازای هر ۱۰ سال 
2022-2102 2 


Climate change with a 0.4°C increase in 
temperature and a 2.15 mm decrease in rainfall 
per decade 
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تغییر اقلیم با کاهش ۴ سانتی گراد دما و افزایش ۲/۱۵ میلی‌متر 


بارندگی به ازای هر ۰ سال 


Climate change with a 0.4°C decrease in 2022-2102 


temperature and a 2.15 mm increase in rainfall 
per decade 


سناریوهای انتخاب شده در مدل برای سال‌های ۲۰۲۲ الی ۲۱۰۲ شامل حفظ وضعیت کنونی (سناریو ۱ 
تغییر اقلیم با افزایش دمای هوا و کاهش میزان بارش به ازای هر ۱۰ سال (سناریو )٢‏ و تغییر اقلیم با کاهش 
دمای هوا و افزایش میزان بارش به ازای هر ۱۰ سال (سناریو (Y‏ هستند. 

Roth C -کالیبراسیون مدل‎ ۴-۳-۲ 

با تغییر فاکتورهای اقلیمی bs)‏ و بارندگی) که از ورودی‌های مدل می‌باشند و نیز نوع پوشش خاک 
ضرایب تجزیه مواد آلی مرتبط با دما و پوشش خاک به‌صورت خودکار درون مدل تغییر می‌يابد. 
بنابراین با وارد کردن اطلاعات اقلیمی مربوط به طبقه ارتفاعی اول در مدلء سرعت تجزیه ماهیانه وابسته 
به دماء رطوبت خاک و وضعیت پوشش گیاهی برای منطقه موردمطالعه تغییر می‌یابد. از آنجا که میزان 
کربن آلی خاک به‌صورت مستقیم اندازه‌گیری می گردد این پارامتر در کالیبراسیون مدل Roth C‏ از 
اھمیت ویژهای برخوردار است. با وارد کردن میزان کربن al‏ خاک در مدل میزان کربن ورودی از طریق 
بقایای گیاهی توسط مدل تخمین و شبیه‌سازی می‌شود. با تغییر میزان کربن ورودی به خاک از طریق SL‏ 
گیاهی, مدل را به‌گونه‌ای کالیبره می‌کنيم که میزان کربن آلی تخمین زده شده توسط مدل با میزان کربن آلی 
اندازه‌گیری شده برابر یا تقریباً مساوی شود. 

Roth C ۵-۳-۲-اعتباررسنجی مدل‎ 

برای اعتبارسنجی» مقادیر ذخیرہ کربن آلی خاک شبیه‌سازی شده توسط مدل با مقادیر محاسبه شده در مناطق 
ارتفاعی دوم سوم و چهارم مقایسه شد و برای انجام این مقایسه و OUS‏ دادن میزان انحراف مقادیر 
شبیه‌سازی شده و محاسبه شده از رگرسیون خطی )2012 (Su, Yan & Tsai,‏ و شاخحص‌های آماری ضریب 
تبیین )2000 (Menard,‏ (رابطه ۳)ء ضریب همبستگی پیرسون )2009 (Cohen et al.,‏ (رابطه (F‏ جذر میانگین 
مربعات (Chai & Draxler, 2014) lke‏ (رابطه ۵) و راندمان مدل‌سازی )2005 (Krause, Boyle & Base,‏ 


(رابطه ۶) استفاده گردید که برای انجام این محاسابات از نرم‌افزار SPSS‏ نسخه ۲۲ استفاده شد. 
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مقدار مقدار مقدار 
ile ; ogy‏ 
á is di Altitude‏ 
Silt Clay Sand (m)‏ 
)96( )96( )96( 
رابطه Y‏ 
رابطه ۴ 
رابطه ۵ 
رابطه ۶ 


۴-۲-آنالیزهای آماری 


آلی n‏ 5 هدایت الکتریکی 
بافت خاک BT Organic‏ اسیدیته Electrical‏ 
or Soil bulk ;‏ 

Conductivity PH densit Carbon Soil Texture 
(ds/m) ۶ ٤ Stocks 
gr/cm (t/ha) 

n ۳ ۳ 2 
R2 = [jt [(6-:0()0-:ک)‎ 


(Oj-0)?‏ و22( 6) بلاط 


(0-:5-0()0) رھ 


۴ VE, (SiS, io] 


n .- ©2 
RMSE =~" x ZrO 


Zi 1(O0i-Si2 


ME =1 - ہج‎ 


بەمنظور بررسی SU‏ ارتفاع منطقه بر ذخیرہ کربن آلی خاک از طرح آماری بلوک‌های کامل تصادفی در سه 


تکرار استفادہ گردید. تجزیه واریانس داده‌ها با استفادہ از نرم‌افزار SPSS‏ نسخه ۲۲ انجام شد )2024 Field,‏ 


۳-نتایج و بحث 


اطلاعات مربوط به pole‏ شن» سیلت و رس. نوع بافت. ذخیرہ کربن آلی؛ چرم مخصوص ظاهری. اسیدیته 


و هدایت الکتریکی برای هر چهار طبقات ارتفاعی در جدول ۶ OUS‏ داده شده است. خاک‌های جنگلی 


موردمطالعه دارای ویژگی‌های یکسانی بودند و از لحاظ میزان رس و نوع بافت. مشابه بودند. همگی فاقد 


آهک بودند و واکنش اسیدی داشتند و مقدار شوری آن‌ها کم بود. 


جدول ۶- اطلاعات مربوط به مقادیر شنء سیلت و رس. نوع بافت» ذخیرہ کربن آلی» چرم مخحصوص ظاهری. اسیدیته 
و هدایت الکتریکی برای هر چهار طبقه ارتفاعی موردمطالعه 


Table 6- Data on sand, silt, clay content, texture type, organic carbon stock, bulk density, 


acidity, and electrical conductivity for each of the four elevation classes under study 
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042 5.49 1.29 44.23 Silty Clay Loam 43.40 3945 17.15 500-1000 
0.36 5.44 1.14 67.20 Silty Clay Loam 49.03 37.90 13.07 1000-1500 
0.26 5.29 1.10 79.47 Silty Clay Loam 47.30 39.01 13.69 1500-2000 
0.19 5.19 1.05 97.46 Silty Clay Loam 49.75 39.35 10.90 2000-2500 


نتایج تجزیه واریانس مربوط به اثر ارتفاع منطقه از سطح دریاء بر ذخیرہ کربن آلی خاک‌ها در جدول ۷ و 
نتایج مقایسه میانگین‌ها در جدول ۸ ارائه شده است. نتایج نشان می‌دهد که اثر عامل ارتفاع از سطح دریا بر 
غلظت و ذخیره کربن آلی SE‏ در سطح یک درصد (۳۰/۰۱< ) معنی‌دار است (جدول ۷) و با افزایش ارتفاع 
از سطح دریاء غلظت و ذخیرہ کربن آلی خاک‌ها نیز افزایش یافتند. به‌طوری‌که بیش‌ترین مقادیر آن‌ها مربوط 
به طبقه ارتفاعی ۲۰۰۰-۲۵۰۰ متری به میزان / ۲۶۲ و ۹۷.۴۶ تن بر هکتار و کم‌ترین daol palis‏ مربوط به 
طبقه ارتفاعی ۵۰۰-۱۰۰۰ متری به میزان / ۰.۹۴ و ۴۴.۲۳ تن بر هکتار بود (جدول ۸). نتیجه این پژوهش با 
نتیجه بسیاری از پژوهش‌های دیگر در مناطق مختلف be‏ هم‌خوانی دارد ) Dieleman, Venter,‏ 
Ramachandra, Krockenberger & Bird, 2013; Banday, Bhardwaj & Pala, 2019; Zhang et al., 1‏ 
یک پژوهش در مغولستان دلیل افزایش کربن آلی خاک در ارتفاعات بالاتر را ناشی از افزایش زیستتودہ 
گیاهی در اثر بارش‌های بیشتر بیان کرده است )1997 (Smith et al.,‏ در جدول ۹ نیز نتایج آزمون همبستگی 
پیرسون بین ارتفاع از سطح دریا و بارندگی منطقه با غلظت و ذخیره کربن آلی خاک ارائه شده است و 
همان‌طور که ملاحظه می‌شود بین غلظت و ذخیره کربن آلی خاک با ارتفاع منطقه از سطح دریا (٭٭0.993- 
ہ و 20.9933 (r‏ و میزان بارندگی ) ٭٭0.958۔ و 0.960% =( همبستگی مثبت و معنی‌داری وجود 
دارد. در پژوهشی دیگری که در ارتفاعات شمال غربی هیمالیا انجام شد. افزايش میزان کربن ST‏ خاک با 
افزایش ارتفاع را به‌سرعت پایین تجزیه مواد آلی در ارتفاعات نسبت داده شد )2019 Banday et al.,‏ به دلیل 
کاهش دما در ارتفاعات بالات فعالیت میکروبی SLS‏ شده و درنتیجه تجزیه مواد آلی و هدر رفت آن کاهش 
می یابد )2007 (Yang, Mohammat, Feng, Zhou & Fang,‏ همان‌طور که در جدول Y‏ نشان oala‏ شده است؛ 
در این پژوهش نیز با افزایش ارتفاع از سطح دریاء متوسط cleo‏ سالیانه منطقه کاهش می‌پابد. به‌طوری که 
بیش‌ترین مقدار متوسط دمای سالیانه مربوط به طبقه ارتفاعی ۱۰۰۰-۵۰۰ متری به میزان ۲۶.۹۳ سانتی‌گراد و 
کم‌ترین مقدار آن مربوط به طبقه ارتفاعی ۲۵۰۰-۲۰۰۰ متری به مقدار ۱۵۸۳ سانتی‌گراد است که این کاهش 


دما در ارتفاعات بالات باعث کاهش سرعت تجزیه مواد آلی و درنتیجه YL‏ رفتن مقدار ذخیرہ آن می‌شود 
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Li, Wang, Endo, Zhao & Kakubari, 2010; Deng, Liu & Shangguan, 2014; Njeru et al., 2017)‏ 
جدول مقایسه میانگین‌ها (جدول ^( OLE‏ می das‏ بیش‌ترین غلظت و ذخیرہ کربن UE‏ در خاک‌های واقع 
شده در طبقه ارتفاعی ۲۰۰۰ تا ۲۵۰۰ متری به میزان ۲۶۲ درصد و ۹۷۴۶ تن بر هکتار اندازه‌گیری شد که 
به‌طور معنی‌داری بیش تر از غلظت و oped‏ کربن آلی اندازه‌گیری شده در ple‏ طبقات ارتفاعی بود. با کاهش 
ارتفاع اراضی از سطح دریا غلظت و ذخیره کربن آلی خاک‌های جنگلی به‌طور معنی‌داری کاهش یافت. 
Sica shes‏ کم‌ترین غلظت و ذخیره کربن آلی به میزان ۰.۹۴ درصد و ۴۴.۲۳ تن بر هکتار در طبقه ارتفاعی 

۰ تا متر اندازه‌گیری شد (حدول MA‏ 

جدول ۷- نتایج تجزیه واریانس SU‏ ارتفاع از سطح دریا منطقه بر غلظت و ذخیره کربن آلی و ذخیره کربن آلی 
خاک‌های جنگلی در منطقه موردمطالعه. 


Table7- Results of variance analysis on the impact of altitude above sea level on organic carbon 
concentration and stock in forest soils in the study area 

میانگین مربعات 

Mean square 


منبع تغییرات 
S.0.V‏ درحه آزادی درصد کربن آلی خاک ذخیرہ کربن الی 
Organic Carbon Organic Carbon Content df‏ 
Stocks‏ 
ار تفا 
e‏ 3 * 1.487 1479.223** 
Altitude‏ 
خطا 
6 0.026 0.209 
Error‏ 


jb‏ 18 به ترتیب معنی‌دار در سطح احتمال پنج درصد. یک درصد و غیرمعنی‌دار. 
جدول ۸- مقایسه میانگین‌های بررسی SU‏ ارتفاع منطقه بر غلظت و ذخیره کربن آلی خاک. 


Table 8- Comparison of means regarding the effect of altitude on organic carbon concentration 


and stock in soil 


ار تفاع 
ویژگی‌های Altitude SE‏ 
Soil properties‏ 500-1000 1000-1500 1500-2000 2000-2500 
صد کرب Í‏ 
درصد کرین الى 0.944 1.68° 2.05 262۰ 
Organic Carbon Content (%)‏ 
"e‏ " 1 
ذخیرہ کربن الى 44.234 67.20° 79.47 ?97.46 


Organic Carbon Stocks (t/ha) 
حروف متفاوت نشان‌دهنده اختلاف معنی‌دار در سطح احتمال پنج درصد در آزمون دانکن است.‎ 


pi ۶۸‏ 4 جغرافیا و مخاطرات محیطی, جلد ۱۳ء شماره ۲ پائیز ۱۴۰۳ 


جدول ۹- ضرایب همبستگی پیرسون (۴) برای ارتفاع از سطح دریاء بارش و دما با غلظت و ذخیرہ کربن آلی 
اندازەگیری شده در موردمطالعه. 


Table 9- Pearson correlation coefficients (r) for altitude, precipitation, and temperature with 
measured organic carbon concentration and stock in the study area 


درصد کربن آلی ذخیرہ کربن آلی 
Organic Carbon Organic Carbon‏ 
Stocks Content‏ 
ارتفا Kok koke‏ 
C‏ 0.993 "0.993 
Altitude‏ 
ہارڈ 
ne‏ "0.958 "0.960 
Precipitation‏ 
m‏ *0.970 *0.971- 
Temperature‏ 


8ء ٭٭ و 18 به ترتیب معنی‌دار در سطح احتمال پنج درصد. یک درصد و غیرمعنی‌دار. 
۱-۳-کالیبراسیون و اعتبارسنجی مدل 
نتایج به‌دست آمده از کالیبراسیون مدل برای ذخیره کربن آلی خاک در اولین طبقه ارتفاعی منطقه موردمطالعه 
در شکل ۳ نشان داده شده است. همان‌طور که مشاهده می‌شود ذخیره کربن آلی شبیه‌سازی شده توسط مدل 
با مقدار اندازه‌گیری شده آن تقریباً برابر است. ذخیره کربن آلی خاک در اولین طبقه ارتفاعی منطقه 
موردمطالعه به میزان» ۴۴.۲۳ تن بر هکتار اندازه‌گیری شده است و مدل نیز مقدار ۴۷۷۳ تن بر هکتار را 


شبیه‌سازی کرده است که مدل مقدار ۳.۵ تن بر هکتار کربن آلی بیشتری را شبیه‌سازی کرده است (شکل (Y‏ 
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شکل ۳- خروجی مدل برای ذخیرہ کربن آلی در طبقه ارتفاعی اول منطقه موردمطالعه تا سال ۲۱۰۲. 


Fig.3. Model output for organic carbon stock in the first elevation class of the study area until 
2102. 


میثاقی و همکاران» شبیه‌سازی تغییرات کربن آلی خاک‌های جنگلی در اثر تغییر اقلیم ... ۶۹ 


ذخیرہ کربن آلی خاک مربوط به سه‌طبقه ارتفاعی دیگر بعد از اعمال کالیبراسیون به مدل داده شد تا دقت 
مدل در شبیه‌سازی ذخیرہ کربن آلی SE‏ مورد ارزیابی قرار گیرد. ذخیرہ کربن آلی (تن بر هکتار) 
اندازەگیری شده برای طبقات ارتفاعی دوم» سوم و چهارم به ترتیب ۶۷.۲۰ ۷۹.۴۷ و ۹۷.۴۶ و مقادیر 
شبیه‌سازی شده توسط مدل برای این طبقات ارتفاعی به ترتیب ۸۷۰.۱۵ ۸۳۵۲ و ۱۰۲.۲۱ است. بدین ترتیب» 
مدل برای طبقات ارتفاعی دوم. سوم و چهارم به ترتیب ۰۲.۹۵ ۴.۰۵ و ۴۷۵ تن بر هکتار کربن آلی بیشتر از 
مقدار اندازه گیری شده. OUS‏ می‌دهد (شکل ۴). 
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شکل ۴- خروجی مدل برای ذخیرہ کربن آلی خاک در هر چهار طبقه ارتفاعی در منطقه موردمطالعه تا سال‎ 
Muy 


Fig.4. Model output for organic carbon stock in soil for all four elevation classes in the study 
area until 2102 


در شکل ۴ء نمودارهای anl‏ 2ص ,m3‏ 104 به ترتیب مربوط به طبقات ارتفاعی اول» دوم سوم و چهارم 
هستند که نشان می‌دهند ذخیره کربن آلی خاک ارتفاعات مختلف با محیط به حالت تعادل است و با گذشت 
ob‏ تغییر نمی کند. 

نتایج مربوط به رگرسیون خطی در شکل ۵ و شاخص‌های آماری ضریب تبیین و ضریب همبستگی پیرسون. 
جذر میانگین مربعات خطا و همچنین راندمان مدل‌سازی که برای اعتبارسنجی مدل لازم می‌باشند. در جدول 


۰ نشان داده شده است. 


pi ۱۷۰‏ 4 جغرافیا و مخاطرات محیطی, جلد ۱۳ء شماره ۲ پائیز ۱۴۰۳ 
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Fig.5. Model validation by comparing simulated and observed organic carbon stock data in all 
four elevation classes in 2022 


جدول ۱۰- شاخصھای آماری محاسبه شده برای اعتبارسنجی مدل. 


Table10- Statistical indices calculated for model validation 


پارامتر شاخص‌های آماری Statistical Indices‏ 
Earameter‏ ضریب همبستگی ضریب on‏ جذر میانگین مربعات خطا راندمان مدل‌سازی 
Modeling RMSE Coefficient of Correlation‏ 
efficiency determination‏ 
کربن آلی 0.97 0.95 4.83 0.94 


Organic Carbon 


نتایج نشان می‌دهد که این مدل می تواند با دقت بسیار بالا ذخیرہ کربن آلی خاک را در منطقه موردمطالعه 
شبیه‌سازی کند. حداکثر اختلاف بین داده‌های ذخیره کربن آلی مشاهده شده و شبیه‌سازی شده ply‏ با ۴۷۵ 
تن بر هکتار (۴.۶۵ درصد) که مربوط به طبقه ارتفاعی چهارم است و کمترین اختلاف نیز برای طبقه ارتفاعی 
دوم به میزان ۲۹۵ تن بر هکتار (۲۰۔۴ درصد) به‌دست آمد. وجود ضریب تبیین (R?)‏ 0.95 نشان می‌دهد ۹۵ 
درصد تغییرات ذخیرہ کربن شبیەسازی شدہ توسط مدل از تغییرات ذخیرہ کربن اندازەگیری شدہ منشأً 
می گیرد و اختلاف واریانس مقادیر اندازەگیری شده و شبیه‌سازی کم می‌باشد )1997 (Smith et al.,‏ ضریب 
همبستگی برای ذخیرہ کربن آلی شبیه‌سازی شده توسط مدل بزرگ‌تر از صفر و برابر با ۰.۹۷ است که 


نشان‌دهنده درجه بالایی از ارتباط و همبستگی مثبت بین polie‏ شییه‌سازی شده و اندازه‌گیری شده است. 


میثاقی و همکاران» شبیه‌سازی تغییرات کربن آلی خاک‌های جنگلی در اثر تغییر اقلیم ... ۱۷۱ 


جذر میانگین مربعات خطا در حد اطمینان ۹۵ درصد برای ذخایر کربن آلی خاک برابر ۴۸۳ بود. محققین از 
جذر میانگین مربعات خطا بەعنوان VEL‏ فاکتور در ارزیابی مدل استفاده کردند ) & Farina, Coleman‏ 
(Whitmore, 2013‏ راندمان مدل‌سازی نیز برابر ۰.۹۴ و مثبت بود که نشان می‌دهد که این مدل دقیق روند 
تغییر ذخایر کربن آلی خاک را در چهار طبقه ارتفاعی موردمطالعه توصیف می کند. 

۲-۳ -شبیه‌سازی سناریوهای مختلف توسط مدل 

سناریوهای انتخاب شده برای مدل برای سال‌های ۲۰۲۲ الی ۲۱۰۲ شامل حفظ وضعیت کنونی (سناریو CV‏ 
تغییر اقلیم با افزايش دمای هوا و کاهش میزان بارش به ازای هر ۱۰ سال (سناریو ۲) و تغییر اقلیم با کاهش 
دمای هوا و افزایش میزان بارش به ازای هر ۱۰ سال (سناریو ۳) می‌باشند. 

۱-۲-۳ -سناریو اول (سناریو بدون تغییر اقلیم) 

در سناریو اول» که حفظ وضعیت کنونی است نتیجه مدل حاکی از در تعادل بودن ذخیره کربن آلی خاک با 
محیط دارد (شکل ۶). به این صورت که مقدار کربن ورودی با خروجی برابر است و ذخیره کربن آلی خاک 
با زمان تغییر نمی‌کند. در شکل ۶ نمودارهای m2 anl‏ 003 و 14 به ترتیب مربوط به طبقات ارتفاعی «Jal‏ 
دوم. سوم و چهارم هستند که نشان می‌دهد ذخیره کربن آلی خاک با محیط در حالت تعادل است و قدرمطلق 
مقدار ذخیره کربن آلی در طبقات ارتفاعی مختلف به دلیل متفاوت بودن شرایط آب‌وهوایی. مختلف است 
ولی در هرسال مقداری دی‌اکسید کربن (CO2)‏ تولید می‌شود که همان کربن خروجی است و مدل OUS‏ 
می‌دهد که مقدار تجمعی دی‌اکسید کربن با گذر زمان افزایش می‌یابد ISS)‏ ۷). علاوه بر این مقدار CO2‏ 
تولیدی در هر طبقه ارتفاعی بسته به شرایط آب‌وهوایی of‏ منطقه متفاوت است و از ۲۷۶.۱ تن بر هکتار در 
سال در بالاترین طبقه ارتفاعی تا ۱۲۸۷ تن بر هکتار در سال در پایین‌ترین طبقه ارتفاعی تغییر می‌کند (شکل 


AV 


۱۷۲ 


نشریه جغرافیا و مخاطرات «dao‏ جلد ۱۳ء شماره ۲ پائیز ۱۴۰۳ 
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شکل ۶- خروجی مدل برای ذخیرہ کربن آلی خاک در هر چھار طبقه ارتفاعی تا سال ۲۱۰۲ (سناریو بدون تغیبر 


اقلیم). 


718.6. Model output for organic carbon stock in soil in all four elevation classes until 2102 (no 


climate change scenario) 
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شکل ۷- خروجی مدل برای دی‌اکسید کربن تجمعی تا سال ۲۱۰۲ در هر چهار طبقه ارتفاعی در منطقه موردمطالعه 


(سناریو بدون تغییر اقلیم). 


Fig.7. Model output for cumulative carbon dioxide until 2102 in all four elevation classes in the 


study area (no climate change scenario) 


میاقی و همکاران. شبیه‌سازی تغییرات کربن آلی خاک‌های جنگلی در اثر تغییر اقلیم ... ۱۷۳ 


۲-۲-۳ -سناریو دوم (سناریو با تغییر اقلیم) 

در سناریو دوم تغییر اقلیم با افزایش ۰.۴ سانتی‌گراد دما و کاهش ۲.۱۵ میلی‌متر بارش به ازای هر ۱۰ سال 
اعمال شد. مدل نشان داد که تا سال ۲۱۰۲ میلادی» میزان ذخیره کربن آلی خاک نسبت به وضعیت کنونی 
(سناریو ۱) کاهش می‌یابد (شکل A‏ بر اساس این سناریو تا پایان سال ۲۱۰۲ میلادی در طبقات ارتفاعی 
اول» دوم. سوم و چهارم به ترتیب ۲ ۴۴.۶۴ ۵۲.۷۶ و ۶۱.۳۳ تن بر هکتار ذخیره کربن آلی برآورد شد 
(شکل (A‏ که به ترتیب ۳۳.۵۴ ۳۶.۳۶ ۳۶,۸۳ و ۴۰۰۰۱ درصد کم‌تر از ذخیره کربن آلی طبقه ارتفاعی مشابه 
در سناریو اول می‌باشد. در شکل A‏ نمودارهای 611 S12‏ 913 و 914 به ترتیب شبیه‌سازی مدل از ذخیره 
ons‏ آلی در طبقات ارتفاعی dal‏ دوم سوم و چهارم با اعمال سناریو دوم است و همان‌طور که ملاحظه 
می‌شود کاهش ذخیره کربن آلی برای LW‏ مرتفع‌تر بیش‌تر است (شکل ۸). نمودار تجمعی دی اکسید کربن 
(CO2)‏ مربوط به سناریو دوم در ٩ JSS‏ نشان داده شده است. نتایج OUS‏ می‌دهد که در سناریو دوم نیز 
مقدار تجمعی دی‌اکسید کربن با گذر زمان افزایش می‌یابد ولی مقدار تولید آن کم‌تر از سناریو اول است که 
این امر به دلیل کاهش میزان کربن ورودی به خاک است (شکل ۹). 
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شکل ۸- خروجی مدل برای ذخیرہ کربن آلی خاک در هر چهار طبقه ارتفاعی تا سال ۲۱۰۲ در منطقه موردمطالعه با 


Fig.8. Model output for organic carbon stock in soil in all four elevation classes until 2102 in 
the study area with climate change (second scenario) 
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Fig.9. Model output for cumulative carbon dioxide until 2102 in all four elevation classes in the 
study area with climate change (second scenario) 


۳-۲-۳-سناریو سوم (سٹاریو با تغییر اقلیم) 

در سناریو cp gan‏ تغیبر اقلیم به‌صورت عکس سناریو دوم اعمال گردید (شکل ۱۰). دال انتخاب این سار 
OUS‏ دادن اهمیت تغییر اقلیم و پارامترهای هواشناسی بر ذخیره کرہن آلی خاک بود. در این سناریو کاهش 
eile‏ کر اه دا و اقرا ۱۵ مل مش رارف بد اڑا هر ١ا‏ مال مایق ھکل ا ر ناس این 
سناریو تا پایان سال ۲۱۰۳ میلادی در طبقات ارتفاعی اول. دوم سوم و چهارم به ترتیب ۵۵.۴۰ ۸۴۳۰ 
۷۰ و ۱۲۶۱۰ تن بر هکتار ذخیرہ کربن آلی برآورد می شود (شکل ۱۰) که به ترتیب ۱۳۸۴ء ۱۶۷۸ء 
۷ و ۱۸.۹۴ درصد بیش‌تر از ذخیره کربن آلی طبقه ارتفاعی مشابه در سناریو اول می‌باشد. در شکل ۱۰ 
نمودارهای 821 822 523 و 824 به ترتیب شبیه‌سازی مدل از ذخیره کربن آلی خاک در طبقات ارتفاعی 
اول. دوم سوم و چهارم با اعمال سناریو سوم است که OUS‏ می‌دهد ذخیره کربن آلی خاک در نقاط مرتفع‌تر 
بیش‌تر افزایش می‌یابد (شکل ۱۰). نمودار تجمعی دی‌اکسید کربن (CO2)‏ تولیدی مربوط به سناریو سوم نیز 


در شکل ۱۱ نشان داده شده است. نتایج نشان می‌دهد که مقدار تجمعی دی‌اکسید کربن با گذر زمان افزایش 


میاقی و ھمکارانء شبیه‌سازی تغییرات کربن آلی خاک‌های جنگلی در اثر تغییر اقلیم ... 


۱۷۵ 


می‌یابد (شکل (VY‏ و مقدار تولید آن نیز بیش‌تر از سناریو اول است که دلیل OF‏ افزایش مقدار کربن ورودی با 
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شکل ۱۰- خروجی مدل برای ذخیرہ کربن آلی خاک در هر چهار طبقه ارتفاعی تا سال ۲۱۰۲ در منطقه موردمطالعه با 


اعمال تغییرات اقلیمی (سناریو سوم). 


Fig.10. Model output for organic carbon stock in soil in all four elevation classes until 2102 in 


the study area with climate change (third scenario) 
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شکل ۱۱- خروجی مدل برای دی اکسید کربن تجمعی تا سال ۲۱۰۲ در هر چهار طبقه ارتفاعی در منطقه موردمطالعه 


با اعمال تغییرات اقلیمی (سناریو سوم). 


Fig.11. Model output for cumulative carbon dioxide until 2102 in all four elevation classes in 


the study area with climate change (third scenario) 


pi ۱۷۶‏ 4 جغرافیا و مخاطرات محیطی, جلد ۱۳ء شماره ۲ پائیز ۱۴۰۳ 


نتایج حاصل از اعمال هر سه سناریو. با یافته‌های سایر محققین مطابقت داشت ) ;2010 Yokozawa et al.,‏ 
ry .(Fallahi, Rezvani Moghaddam, Nassiri Mahallati & Behdani, 2013; Soleimani et al., 2017‏ وھی L‏ 
استفاده از این مدل کاهش محتوای کربن آلی خاک را در سناریوهای مختلف تغییر اقلیم نسبت به سناریو 
بدون تغییر اقلیم برای دهه‌های آیندہ در سراسر اسپانیا پیش‌بینی نمودند ) Jebari, Del Prado, Pardo,‏ 
(Rodriguez Martin & Alvaro-Fuentes, 2018‏ تحقیقات نشان داده است که در Ole‏ پارامترهای اقلیمی» دما 
و بارندگی بیشترین P‏ را ہر مقدار مادہ آلی خاک دارد )2003 (Shakiba & Rahnama,‏ شواهد حاکی از آن 
است که مقدار کربن آلی خاک با کاهش میزان بارندگی کاهش می‌یابد» ولی در هر سطحی از بارش مقدار 
ماده آلی خاک با افزايش clo‏ هوا کاهش می (Singh, Pandey, Sidhu, Sarkar & Sagar, 2011) Jub‏ نتایج 
اکثر مطالعات انجام شده OUS‏ می‌دهد افزايش دمای هوا با افزایش سرعت تجزیه میکروبی و فرایند معدنی 
شدن باعث کاهش مقدار ماده a‏ خاک می‌شود ) Smith et al., 1997; Kóchy, Don, van der Molen‏ 
(Freibauer, 5‏ و تعداد کمتری از محققین نیز بر این باورند که افزایش دمای هوا موجب افزایش میزان 
تولید بقایای گیاهی و درنتیجه باعث افزایش کربن ورودی به خاک و درنھایت میزان ماده آلی خاک می‌شود 
«Rampazzo Todorovic et al., 2010)‏ که با یافته‌های این پژوهش هم‌خوانی ندارد. 

در شکل ۱۲ به‌طور خلاصه مقدار کربن آلی ذخیرہ شده در هر چهار طبقه ارتفاعی در هر سه سناریو اعمال 
شده نشان داده شده است و این شکل OUS‏ می‌دهد در تمام سناریوها نقاط مرتفع‌تر» دارای ذخیره کربن آلی 
بیشتری نسبت به نقاط پست‌تر هستند. علاوه بر این در اثر تغییر اقلیمی که با کاهش بارندگی و افزایش دما 
همراه می‌باشد ذخیره کربن آلی خاک‌های جنگلی طبقات ارتفاعی مختلف کاهش می‌یابد. درصورتی که 
چنانچه تغییر اقلیم با کاهش دما و افزایش بارندگی همراه باشد ذخیره کربن آلی خاک افزایش می‌یابد و این 


افزایش برای نقاط مرتفع‌تر بیشتر است. 


میاقی و ھمکارانء شبیه‌سازی تغییرات کربن آلی خاک‌های جنگلی در اثر تغییر اقلیم ... 


Organic carbon stock in soil under all three scenarios 
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Fig.12. Amount of organic carbon stock in each of the four elevation classes in the study area 


for all three applied scenarios. 


۴-نتیجه گیری 


در این مطالعه تأثیر ارتفاع از سطح دریا و تغییر اقلیم بر ذخیرہ کربن آلی خاک‌های جنگلی در منطقه تالش 


استان OAS‏ موردبررسی قرار گرفت. نتایج نشان داد که با افزایش ارتفاع از سطح دریاء ذخیرہ کربن آلی 


خاک‌های جنگلی افزایش می‌یابد. این افزایش به دلیل عواملی مانند افزایش میزان بارندگی و کاهش دما 
است. تغییر اقلیم نیز می‌تواند بر ذخیرہ کربن آلی خاک‌های جنگلی تأثیر قابل‌توجهی داشته باشد. 


درصورتی که تغییر اقلیم منجر به افزايش دما و کاهش بارندگی شود ذخیره کربن آلی خاک‌های جنگلی 


کاهش خواهد یافت. در مقابلء اگر تغییر اقلیم منجر به کاهش دما و افزایش بارندگی شود. ذخیره کربن آلی 


خاک‌های جنگلی افزايش خواهد Cob‏ این مطالعه نشان داد حاک‌های جنگلی, به‌ویژه خاک‌های واقع شده 


در ارتفاعات SVL‏ نقش مهمی در جذب و ذخیرہ ous‏ اتمسفری دارند. تغییر اقلیم می تواند این نقش را 


تحت تأثیر قرار دهد؛ بنابراین. مدیریت پایدار جنگل‌ها و خاک‌های جنگلی» برای کاهش اثرات تغییر اقلیم بر 


محبط زیست. ضروری cual‏ در مطالعات آینده» AL‏ به کمک مدل‌های مو جود. تأثیر اقدامات مدیریتی که 


باعث کاھش سرعت تجزیه ماده آلی و هدر رفت آن به‌صورت گاز کربنیک می شود موردبررسی قرار گیرند. 


۱۴۰۳ جغرافیا و مخاطرات محیطی, جلد ۱۳ء شماره ۲ پائیز‎ 4 pi VA 


همچنین» در دسترس نبودن داده‌های هواشناسی برای خاک‌های جنگلی واقع در ارتفاعات» یکی از 


محدودیت‌های این پژوهش است که در آینده باید به آن توجه شود. 
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